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I.  INTRODUCTION 


The  whole-body  motion  of  an  item  of  military  equipment  struck  by  a 
long-duration  blast  wave  may  result  in  functional  or  physical  damage  to 
that  item.  In  order  to  develop  techniques  for  predicting  such  motion,  it 
is  necessary  to  describe  the  loading  on  the  target  produced  by  the  blast 
wave.  In  a computer  code  developed  at  the  Ballistic  Research  Laboratory 
to  calculate  rigid  body  rotation,  the  target  was  modeled  as  a collection 
of  rectangular  parallelepipeds^.  For  such  a code,  a load  prediction 
technique  is  required  which  considers  the  variations  in  dimensions  of 
such  parallelepipeds  which  occur  in  a given  model. 

The  blast  loading  of  a target  occurs  in  two  phases;  the  diffraction 
phase,  where  the  shock  wave  encounters  and  engulfs  the  target;  and  the 
drag  phase,  which  occur.s  when  variations  due  to  the  diffraction  process 
h.ive  ceased  and  quas i -steady  drag  loading  is  experienced  by  the  target. 

This  report  is  concerned  with  developing  a prediction  technique  for  the 
average  pressures  developed  on  the  front  and  rear  surfaces  of  a rectang- 
ular parallelepiped  during  the  diffraction  ph.nse.  The  target  is  assumed 
to  bo  struck  by  a shock  front  which  is  moving  nomal  to  the  front  surface. 

No  satisfactory  diffraction  loading  prediction  technique  was  available 
for  predicting  the  average  loading  on  all  surfaces  of  a target.  In  Ref- 
erence 2,  Taylor  describes  a comparison  between  results  obtained  in  shock 
tube  ex{)criments  and  loading  predicted  by  techniques  contained  in  .standard 
manuals''.  Significant  differences  wore  shown,  and  therefore  such  techni- 
ques require  iapi'oveuent . 

The  basic  data  used  by  Taylor  for  the  analysis  described  in  Reference 
2 wore  available  at  BRl** . The  purjwse  of  this  report  is  to  present  equa- 
tions fitted  to  these  data  to  provide  a prediction  technique  for  the  aver- 
age loading  on  the  front  and  rear  t.^rfaces  of  rectangular  parallelepipeds. 
l:qu,atlons  arc  prcsentisl  as  fitted  to  the  data  and  also  as  liiodified  for 
use  in  a rigid-body  response  code. 


^*,V.  H.  Sthvu^Cf  "BLiat  fJiwrSurntnj?  /or  farvteto,"  BRl 

AVfKJi’t.  im.  (AD  iB0i2“^l02L) 

"V,  */.  Taylor,  "d  Hgthod  for  PiSidiotiHrs  flZtsne  boadu  thtritiy  tho  sHffmoiion 
Phiso, " USA  Rdocatsih  fjdboiMtori^sa,  Abanieen  Pravi*^y  Ut'ound, 

^tai'yland;  inaludad  in  Birt  4,  Bulletin  45,  ffce  Sltaok  ami  Vibratim 
Balltitin,  The  Shack  and  Vibrntian  Infomstim  Center,  liaml  Researoh 
Laboratory,  ihahifujUm,  0.  C.,  January  i$75. 

t 

C.  H.  Rorrio,  et  al,,  **Strmtural  Ueoign  for  Dymmio  Loada,*'  HeGimt-liiti 
lk>ok  Con^my,  Ino.,  Neu  York,  B.  Y,,  1559. 

Taylor,  Prioate  Comutiioatim. 
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Figure  1 shows  the  target  used  by  Taylor  in  the  shock  tube  experiments. 
Pressure  gages  were  positioned  at  locations  labeled  A,  C,  and  D.  Position 
B was  not  instrumented,  since  by  symmetry  the  pressure  at  Position  B should 
be  the  same  as  that  at  Position  A.  To  obtain  the  average  pressure  the  out- 
put of  all  three  gages  were  electronically  summed.  The  output  from  the 
gage  at  A was  given  twice  the  weight  of  the  output  of  gages  D and  C to 
effectively  include  the  record  that  would  have  been  obtained  by  a gage  at 
B. 


Figure  2 shows  the  test  configurations  used  by  Taylor.  By  the  addition 
of  blocks  on  the  ends  of  the  model  a two-dimensional  configuration  was  pro- 
duced. Data  were  obtained  for  both  configurations  and  for  front  and  rear 
surfaces  for  nominal  shock  overpressures  of  54.5,  69,  105,  and  158  kPa  (5, 
10,  15,  and  20  psi).  Tables  i and  2 list  the  tost  shots  and  data  appro- 
priate for  each^ . 

Gentry^  made  a hydrocode  calculation  for  this  target  for  a shock  wave 
of  54.47  kPa  (5  psi)  incident  overpressure  for  the  three-dimensional  con- 
figuration shown  in  Figures  2(A)  and  2{C).  The  hydrocode  predictions  for 
individual  gage  records  matched  the  experimental  data  very  well.  Average 
pre.ssures  computed  for  the  front  and  roar  surfaces  wore  in  good  agreement 
with  average  pressures  calculated  using  computed  values  at  the  gage  posi- 
tions and  the  same  averaging  procedure  as  was  used  for  the  gage  records 
in  Reference  2.  Thu.s  the  procedure  used  by  Taylor  to  record  average  pres- 
sure on  front  and  rear  surfaces  sotms  satisfactory  for  the  three-dimen- 
sional configuration.  No  hydrocode  calculation  was  available  to  check  the 
averaging  preveedure  for  the  two -dimensional  configuration  shown  in  Fig- 
ures 2(B)  and  2(1)).  The  positions  of  the  gages  and  the  averaging  proce- 
vlurc  do  not  sgwa  as  appropriate  for  the  two-dimensional  configuration  as 
for  the  three-dimonslonaS  configuration. 


11.  FRONT  SUUFACK  ld3iAi)I.NG 

A.  Fitting  an  Equation  for  Average  Front  Surface  Pressure 

In  Reference  2,  Tsylor  found  that  two  and  three-dimensional  average 
pressure  curves  on  the  front  surface  of  the  shock  tube  target  would 
nearly  coincide  if  the  time  scale  were  described  in  units  of  the  number 
of  rarefaction  wave  crossings  of  the  loaded  area.  A rarefaction  wave 
originates  at  each  free  edge,  so  that  for  the  three-dimensional  con- 
figuration two  rarefaction  wave  crossings  occur  in  the  time  that  one 
occurs  for  the  two-dimensional  configuration.  Figures  .5  and  4 .show 
figures  from  Reference  2 which  present  the  average  front-surface  pressure 
loading  versus  the  number  of  rarefaction  wave  crossings.  These  data  were 
used  for  fitting  an  equation.  M noted  above,  the  thrce-diacnsional  data 


F,  A.  (ktKtiVf,  iii.,  Dimn^ional  CoKputBr  Arushjaia  of  Skottk 

Iciido  on  n Sv^^le  StJ'itatum,”  BHl  CF  S25,  VSA  Balliatia 
Ltabomtoin^io,  Frooing  Gvound^  fSaiSih  1$?6.  (Al)  #80032081) 
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(A)  SURFACE  ZONING 


Figure  1 Target  diinenalons  and  gage  positioning  used 
in  the  study  described  in  Reference  2. 
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apparently  represent  the  average  pressure  very  well,  while  the  two-dimen- 
sional data  may  not  present  an  accurate  average.  However,  in  the  absence 
of  other  data  the  two  dimensional  curves  were  fitted  and  used  in  forming 
a prediction  equation. 

The  basic  form  of  the  equation  chosen  for  the  fitting  was  as  follows; 


(1) 


where : 

P = average  overpressure 

P|^  = peak  reflected  overpressure 

P^  = overpressure  approached  as  a limit 

T St  a^Cl  ♦ R)/S 

t = time 

85  » velocity  of  the  rarefaction  wave  behind  the  reflected  shock 

wave 

R e ratio  of  the  odgc.s  of  the  minimum  rectangular  area  which  must 
be  examined  to  define  the  loading  over  the  entire  surface  area. 
The  ratio  is  taken  so  that  R ranges  from  0 for  tlie  two-dimen- 
sional configuration  to  1 for  the  three-dimensional  configura- 
tion. 

S « minimum  distance  of  travel  from  an  edge  for  a rarefaction  wave 
to  traverse  the  loaded  area. 

A,  B a functions  of  time,  incident  ovei'pressuro,  ambient  pressure, 
and  R. 

This  equation  provides  the  appi’opriate  tero  Initial  slope  and  can  be 
made  asymptotic  to  the  stagnation  overpressure  or  some  other  selected 
pressure  limit. 

The  limiting  average  overj>ressure  for  P,  was  chosen  to  be  the  stagna- 

g ^ 

tion  overi>ressurc,  Actually  the  average  limiting  overpressure  on 

the  front  surface  should  be  somewhat  less  than  because  some  flow 

* occurs  across  the  surface  as  the  air  flows  around  the  target.  In  Hefet- 

enco  S,  the  average  overpressure  limit  is  given  as  P + 0.8S  P , where 

^ H 

P^  is  the  incid«uit  peak  overpressure  and  P^  is  the  correspmidittg  dynamic 

pressure.  One  would  expect  the  aiwnint  of  the  reduction  from  stagnation 
overpressure  to  be  a function  of  the  shape  parameter  B,  P^,  and  P^.  i 
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The  limited  set  of  data  available  from  Reference  2 suggest  that  the 
three-dimensional  data  may  indeed  approach  a limit  less  than  that 
approached  by  the  two-dimensional  data.  However,  because  of  the  limited 
data,  was  chosen  as  the  average  overpressure  approached  as  a limit 

for  this  fitting  exercise. 

The  function  A was  found  by  examining  values  read  from  the  curves  for 
T = 1.  For  that  value  the  equation  becomes: 


Solving  this  equation  for  A,  the  result  is: 


Data  for  the  two  and  three-dimensional  configurations  at  t = 1 were 
found  to  bo  almost  equal  for  the  same  incident  shock  ovorpressuro , so 
those  data  wore  combined  and  fitted  as  a function  of  z,  where  z is  the 
ratio  of  incident  shock  overpressure  Fg  to  the  ambient  atmospheric  pres- 
sure I*  . 

A 

The  function  found  for  A was  inserted  into  Equation  1 and  corresponding 
values  of  the  fvmction  B wore  calculated  to  provide  a smooth  representation 
of  the  curves  in  Figures  .>  and  4.  A fimction  was  fitted  to  those  values 
of  B versus  T.  A constant  value  provided  a satisfactory  reprosontaiton  of 
the  three-dimensional  data.  However,  for  the  two-dimensional  data  a vari- 
ation with  both  2 and  t was  found  noces.sary  for  a reasonable  representation 
and  th.  form  selected  was: 

B iz,  T)  « a ♦ e ^ (S) 

where  a,  a and  8 were  the  chosen  parameters, 

Figure  5 shows  the  area  on  the  front  surface  of  a block  target  for 
which  pressure  loading  must  be  determined  in  order  to  define  the  loading 
over  the  entire  front  surface  area.  If  the  target  is  in  free  air  there 
are  two  planes  of  symmetry,  and  only  one-fourth  the  area  need  be  examined. 
For  a target  on  the  surface,  which  represents  the  configurations  shown  in 
Figures  2A  and  2C,  a vertical  plane  of  sywsetry  exists  and  only  one-half 
the  area  must  examined.  The  parameter  R is  defined  as  the  ratio  of  the 
two  edges  of  the  ttinimus  rectangular  area  which  must  be  exainin<^  to  diefinc 
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on  aurface  agalnat 
vertical  wall 


on  aurface 


Area  used  for 
CoapuUtion 

in  free  air 


Fiaure  5 Minimum  rectatiRUlar  area  whlcb  nust  be  exaalucd  to 
ngur.  5 ^ ^ surf.™  for  sovera 

oxpoQure  conditions • 


pressure  loading  over  the  entire  front  surface.  The  ratio  is  taken  so 
that  R is  less  than  or  equal  to  one.  For  the  two-dimensional  data  ob- 
tained for  the  configurations  shown  in  Figures  2B  and  2D  the  value  of  R 
is  zero.  For  the  three-dimensional  configurations  shown  in  Figures  2A 
and  2C  the  value  of  R is  essentially  one. 


It  was  assumed  that  the  average  pressure  loading  would  vary  from  that 
fcr  R = 0 to  that  for  R = 1 in  proportion  to  the  value  of  R.  In  the 
fitted  equations,  then,  the  parameter  R is  inserted  whore  a difference 
due  to  configuration  occurs  to  produce  a smooth  variation  with  the  value 
of  R.  The  final  form  of  the  fitted  equation  for  front-surface  loading 
was  'Equation  1 with; 


where : 


A = 
B B 


P 

stag 


0.265  z-0.20 

„ ,0.37 

1.5R  ♦ (1  - R)  {0.97  ♦ ^ (4) 

2 Pg(4z  ♦ 7)/(z  + 7) 

P . 

stag 

stagnation  overpressure  on  the  front  surface;  which  for 
Mach  number  less  than  ono  is: 


P . = (P  ♦ P jfl  ♦ 

stag  s A'l  7(2+1;{s*7)J 


P = incident  shock  overi)rossure 


P^  “ ambient  .atmospheric  pre.ssuro 


z - P./P^ 


7/2 


*>  sound  velocity 


- p. 


a/s 


a ambient  atmospheric  temperature,  degrt'^^s  Kelvin 

Figures  6 through  13  show  the  curves  from  Figures  3 and  4 cDttpareU 
with  values  calculated  using  Equation  1 with  functions  defined  as  shown 
in  Equation  4. 
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1.1 


Figure  6 Average  tuceaaure  ratio  on  the  front  face  of  tlie  two- 
dlmenalonal  shock  tube  taiget  eotapared  with  fitted 
eq^uatlona  and  a statidard  pr6dlctl(m  techolciue  for  an 
incident  shock  overticressure  of  35*13  (5*24  pal). 
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p/p. 


72.2  kl^  (10. U?  psi) 


r 

Figure  7 Average  pressure  ratio  aa  the  front  face  of  the  two- 
dlaeneional  shock  tube  target  ccaapared  with  fitted 
eiiuatiojiB  and  a staiidard  prediction  technique  for  an 
incident  shock  ovcrpreasure  of  7^.2  kBa  (10.47  psi). 
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T 

Pleuife  8 Avejrage  pceasure  ratio  oa  the  front,  face  of  the  two- 
dimenalc«al  uhocii  tube  target  cmparcd  with  fitted 
equatloaa  end  a standard  isredletion  technique  for  an 
incident  chock  overi»reiicure  of  108^9  kBa  (15.3  psi). 


- ' 


34.5  kPa  (5.0  psi) 

— Sh-.'t  1 - S-DiniRi’Sional 

— Equations  1 and  4 

-X  Standard  Prediction  Technique 

r = t *•  B)/^ 
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Ilgure  10  Average  prosaure  ratio  on  th:'  front  face  of  the  tl>ree- 
dimenslonal  shock  tube  target  compai’t<?  with  fitted 
equations  and  a standard  preductlon  teclmique  for  an 
incident  shock  cverpreasiu’e  of  3*^»5  hPa  (5.0  psi). 
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dimensional  shock  tube  target  ccsapared  with  fitted 
equations  ajid  a etandaM  prediction  tsehnique  for  an 
Incident  oliock  overpreaaure  of  73,7  k^  (XO.7  pal). 
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Figure  12  Average  treasure  ratio  oa  tiie  front  face  of  the  three- 
diKo^^ional  ohoek  tube  tar^t  ooapared  vith  fitted 
equations  and  a otanilai'd  prediction  technique  for  an 
incident  {fhock  overpreasure  of  1X1.1  kBa  (16.1  psi). 

2& 


^3  Average  preaoure  ratio  on  the  ivont  face  of  the  three 
dimeE^ional  shock  tube  target  compared  with  fitted 
e<iuationa  and  a etandard  prediction  technique  for  an 
incident  ahook  ovorpreaaure  of  1^7 •X  kBa  (2X.3  pai)< 


The  quality  of  fit  to  the  experimental  data  was  judged  by  inspection 
rather  than  on  statistical  measures  which  would  be  provided  by  a least- 
squares  fit  of  Equation  1 to  the  data.  The  fits  developed  are  estimated 
to  be  within  the  range  of  error  for  gage  records  and  the  gage  record 
averaging  procedure. 

Before  performing  the  fits,  the  curves  from  Figures  3 and  4 were 
normalized  so  that  the  peak  value  corresponded  to  the  peak  reflected 
pressure.  These  experimental  summed  gage  records  should  show  an  initial 
period  of  no  decay  which  corresponds  to  the  time  for  a rarefaction  wave 
to  reach  the  nearest  gage  from  an  edge.  The  actual  average  loading  curve 
will  begin  to  decay  immediately. 


The  equation  was  fitted  to  data  for  overpressures  ranging  from  34.5  ;1 
to  138  kPa  (5  to  20  psi) . The  equation  can  be  used  outside  this  range,  ■; 
say  from  13.8  to  345  kPa  (2  to  50  psi),  but  with  less  confidence.  i; 
Figure  14  shows  curves  calculated  using  Equation  1 with  functions  as  | 
defined  in  Equation  4 for  13.8,  34.47,  137.9,  and  344.7  kPa  (2,  5,  20,  > 
and  50  psi)  for  = 101.325  kPa  (14.696  psi)  and  T^  = 15°C  and  R = 0.5,  | 

and  Y = 1.4.  I 


B . Comparisons  with  Other  Shock  Tube  Data  j 

Reference  6 presents  shock  tube  data  on  front- surface  loading  obtained  - 

on  a target  with  R » 1.0,  corresponding  to  the  three-dimensional  config-  | 

uration,  for  overpressures  of  22,  40,  62,  and  100  kPa  (3.2,  5.8,  9.0,  ■ 

and  14.5  psi)  at  an  average  P.  of  102  kPa  (14.8  psi).  The  target  was  f 

A .»! 

SO. 8 mm  (2  in.)  high,  101.6  mm  (4  in.)  wide  and  50.8  mm  (2  in.)  deep.  j 

Five  gages  were  located  on  the  front  face.  The  individual  gages  wore  | 

recorded  separately.  The  records  wore  road  and  the  values  combined  to  | 

form  the  average  front-surface  pressure.  The  records  show  a finite  rise  1 

time,  and  the  peak  values  do  not  correspond  to  the  theoretical  values  I 

of  peak  reflected  pressures  appropriate  for  the  values  quoted  for  the  | 

incident  ovorpresauro.  However,  it  is  not  clear  that  accuracy  of  the  | 

average  pressure  at  later  times  would  bo  improved  by  normalizing  the  1 

rocoi'ds  so  that  the  peak  values  shown  in  Figui'os  21-24  of  Reference  6 t 

would  corresjKjnd  to  the  theoretical  peak  reflected  pressure.^.  Those  | 

curves  wore  transformed  to  plots  of  P /P,,  versus  t without  attempt-  I 

average  R >; 

ing  any  normalization.  | 

I 

The  transformed  curves  are  shown  in  Figures  1$  and  16.  The  dashed  | 

curves  wore  calculated  using  Equation  1 with  functions  as  defined  in 
Equation  4.  The  agreement  between  the  predicted  curves  and  the  data 
curves  is  best  for  the  lowest  incident  shock  overpressure.  Agreement 
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13.79  ItPa  (2  psi) 
34.47  kPa  (5  psi) 

|137.9  kPa  (20  psi) 

j 344.7  l£Pa  (50  pel) 


Figure  l4  Average  precsui’e  I'atioa  calculated  usiug  Equations  1 
and  4 for  13.79,  137. 9,  344.7  kia  (S,  5> 

20,  aiid  50  pai)  for  « 101.325  kBa  (14.696  pal) 

aM  T,  « 283.15  de^eea  Kelvin  and  B 0.5. 
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Figure  15  Cottpaapiaon  of  Equations  1 and  4 with  average  jjreesure 
ratios*  <«  the  front  face  of  shock  tube  target  of 
Reference  6.  Incident  shock  overpressures  v<ere  22 
end  1(0  kFa  (3.2  and  5<d  ipai)  and  it  « 1.0. 
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16  Comparison  of  EquatiooB  1 ani  4 vlth  avemfio  l^esure 
ratios  CO  the  front  face  of  the  shocS;  tube  target  of 
Heterence  6,  Incideat  shock  overpressures  were  62 
and  ^00  kBa  <9^0  and  X4.§  psi)  and  B « i.O. 


would  indeed  be  improved  if  normalization  of  the  shock  tube  data  to  the 
theoretical  values  of  peak  reflected  pressure  were  justified.  The  general 
shape  of  the  experimental  curves  seems  to  correspond  to  the  predicted 
curves.  The  data  indicate  a decay  to  an  overpressure  below  the  stagna- 
tion pressure. 

C.  Comparison  with  a Three-Dimensional  Hydrocode  Calculation 

In  Reference  7,  Lottero  presents  the  results  of  a three-dimensional 
hydrocode  calculation  of  the  loading  developed  on  an  S-280  Electrical 
Equiimient  Shelter  by  a 34.47S  kPa  (5  psi)  shock  wave.  The  shock  wave 
was  assumed  to  have  no  decay  in  pressure  behind  the  shock  front,  and  to 
strike  the  side  of  the  shelter  at  normal  incidence.  The  calculation 
was  made  by  R.  A.  Gentry,  et  al.,  at  the  Los  Alamos  Scientific  Laboratory 
using  the  hydrocode  BAAL.  The  calculation  was  made  in  the  same  fashion 
as  that  reported  by  Ctcntry  in  Reference  5. 


Figure  17  shows  a plot  of  predictions  made  using  Equations  I and  4 

compared  with  the  average  pressure  computed  from  the  hydrocode  results. 

The  ratio  of  P /P„  is  shown  plotted  versus  time.  Values  of  the 

average  R 

parameters  used  in  Equation  4 were  ® 0.4904,  S » 1.81  metros, 

R o 0.8584,  a_  = 369.62  m/s,  and  z = 0.3402. 

d 

The  hydrocode  data  indicate  an  initial  decay  to  a pressure  value 
less  than  the  stagnation  pressiu'o  limit  of  Equation  I with  functions 
as  dcfinotl  in  Equation  4.  A pressure  limit  somewhat  loss  than  stagnation 
prossui‘0  is  physically  rea.sonablc,  since  there  will  be  .some  flow  across 
the  front-surface  area  UTid  near  the  edges  the  pressure  will  not  bo  at 
the  full  .stagnation  value.  The  agreement  between  the  pi’odictod  curve 
and  the  hydrocode  data  points  seems  satisfactory. 


0 . Comparison  with  Standard  Prediction  Techniques 


In  Reference  .3  the  average  pressure  on  the  front  surface  of  a target 
is  presented  as  decaying  from  the  peak  reflected  pressure  to  a value 
P^  ♦ 0,85  where  P^  is  the  overpressure  and  P^  is  the  dynamic  pressure 

in  the  blast  wave.  The  decay  is  approximated  by  a .straight  line  and 
occurs  in  the  time  t,  « 3 S/a,^.  This  time  corresponds  to  t “ 3 for 


the  two-dimensional  configuration  where  R » 0 and  to  t 6 for  the 
three-dimensional  configuration  where  R = 1.  In  Figures  6 through  13 
the  predictions  made  using  this  technique  are  shown  by  the  straight  linos 
terminating  at  an  X.  The  disagreement  is  particularly  largo  for  the 
three -dimensional  configuration.  The  primary  purpose  of  the  work  by 
Taylor^  was  to  point  out  this  disagreemont  and  to  emphasize  that  the 
standard  prediction  techniques  were  not  very  good  when  applied  to  the 
threc-dlnensional  configuration. 


B.  S,  '^Coifputatiaml  PnHsdiatione  of  S^took  Diffisiotion  Loading 

f>n  an  S^MO  Bl&otHoal  Sh^ltePt"  BHL  UmosKmdm  Baport  $S09t 

ms.  (AD  »A02a804) 
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Figure  17  OoBjparisoa  of  average  preaeure  ratio  coaiputed^  for  an 
8>£!60  Blectricad  Equlp»nt  Siuilter  xiaing  a tb^ee- 
dloeasioael  hydrocode  (BAAL)  vith  Equations  1 and  4 
f or  au  incident  ehoclt  overpreaaure  of  34.4?  isBa  pal). 


E.  Modification  to  Account  for  the  Decay  of  the  Incident  Blast  Wave 

The  data  obtained  by  Taylor^  apply  for  a shock  wave  with  no  pressure 
decay  behind  the  shock  front,  which  corresponds  to  a blast  wave  from  an 
infinitely  large  explosion.  To  account  for  the  pressure  decay  which 
occurs  in  a real  blast  wave,  the  pressures  predicted  by  Equation  1 are 
reduced  in  proportion  to  the  reduction  in  with  time.  This  provides 

only  an  approximation  to  the  correct  pressure,  since  the  reflected  pres- 
sure is  a non-linear  function  of  the  incident  shock  overpressure  waveform. 
However,  the  correct  function  is  not  available.  A treatment  of  non-linear 
shock  wave  reflection  theory  in  one-dimensional  form  is  given  in  Reference 
8,  and  is  quite  complex,  even  for  one-dimensional  flow. 

If  this  procedure  is  incorporated  into  Equation  1 to  account  for  the 
decay  with  time  of  the  incident  shock  overpressure,  the  prediction  equa- 
tion becomes  as  follows; 

where:  P(t)  » average  front-surface  overpressure  at  time  t 

(^)  “ stagnation  overpressure  in  the  blast  wave  at  time  t 

“ “’s")  ‘ ’’a'  * 7(p;cn-v-iyJ  - '’a 

Pg(t)  ® overpressure  in  the  blast  wave  at  time  t 

P « stagnation  overpressure  corresponding  to  peak  incident 
overpressure  1*^. 

F.  An  Alternative  Pit  to  the  Average  Front -Surface  Overpressure 

In  the  course  of  using  Equation  I for  predictions  it  was  noted  that 

predicted  curves  for  the  extreme  ends  of  the  overpressure  range  of 

interest  (13  to  34S  kPa,  or  2 to  SO  psi)  did  not  differ  very  ouch  in 

initial  shape  when  plotted  as  P ^ „ /P„  versus  ti«e.  This  result 
* average  R 

suggested  examination  of  the  data  in  terms  of  a time  normalisation  with 
no  dependence  on  ovorpressuro , such  os  «•  t a^(l  ♦ R)/S,  whore  is 

the  sbimd  velocity  in  the  undistiurbed  air* 


u 
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Figures  18  and  19  show  the  two  and  three-dimensional  data  for 
^average'^^R  versus  t*.  The  correspondence  is  such  that  no 

attempt  was  made  to  maintain  the  separate  identity  of  the  curves  in 
the  initial  part  of  the  plots.  There  is  no  indication  of  a consistent 
variation  with  overpressure  level  until  most  of  the  pressure  decay  has 
occurred. 


Equations  of  the  same  basic  form  as  Equation  1 were  fitted  to  the 
data.  For  the  two-dimensional  case,  R = 0,  the  results  for  the  functions 
A and  B were: 

A = 0.388 

B = l.lSz'^’^  + 0.026z^‘^^  (t-6)^  i6) 


where  x = t*  » t ag  (1  ♦ R)/S,  and  other  quantities  have  been  defined  for 
Equations  1 and  4. 

For  the  three-dimensional  data  the  fit  was  made  neglecting  the  appar- 
ent variation  in  the  limiting  pressure  at  late  times.  The  stagnation 
pressure  was  used  as  the  limit  to  simplify  joining  the  resulting  equation 
with  Equation  6 to  obtain  one  equation  applying  for  a'.l  values  of  P.  between 
0 and  1.  The  results  for  A atwl  D for  the  three-dimensional  configuration 
were : 

A s 0.M(3 

8 = 1.54  ♦ 0.1  (x-2.3)”  (7) 


P.  a P 
L stag 


T “ X* 

The  final  functions  combining  the  results  of  both  fittings  are: 

A » 0.388  (i-R)  ♦ 0.346  R 

8 « (1-R)  [i.l3t"®*^  ♦ 0.0262®*^^  {t-6)“]  + R [1.54  ♦ 0.1  (t-2,5)^3  (8) 


t « X* 

Figures  20  through  37  show  the  curves  calculated  from  Equation  1 with 
functions  as  defined  in  Equation  9 compared  with  the  experimental  curves. 
It  is  clear  that  the  resulting  curves  represent  reasonable  fits.  Equa- 
tion 1 can  be  used  with  functions  as  defined  by  either  Equation  4 or 
Equation  8 in  the  range  from  34. S to  138  kPa  (S  to  20  psi)  with  no 
ficant  difference  in  the  result. 
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PlguM  iB  Avex-age  ftront  surface  pressure  ratios  of  the  two- 
diiiiettSlonal  ahoefc  tube  target  of  Befereuce  2 
plotted  versus  ^ * ^)/S’ 


19  front  surfnse  prcsaikre  ratios  of  tiie  three 

dlsasc^lpoiU.  ahock  tube  target  of  Kefereiice  2 
plotted  versus  t a_(l  ♦ R)/S, 


1.1 


0.3  ir-'i- 

I ■ 
f-'-i;  , 
ir.Ns;,;  ^ 


36.13  kPa  (5.24  psi) 

— Shot  17  - 2-J)imen5ional 

— Equations  1 and  8 

n = t a (1  + R)/S 


If  .t.  Vi 
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Figure  20  Average  front  surface  pressvire  ratio  of  the  t'^o-dlnranalonal 
shock  tube  target  for  aft  incident  shock  overpressure  of 
36.13  kPa  (5.24  psi)  plotted  vers’is  ^ » t a {l  + K)/S. 


Figure  22  Average  front  surface  pressure  ratio  of  the  two-dimensional 
shock  tube  teu'got  for  m incident  shock  overpwaaure  of 
108.9  Itfia  Ci5.d  p«i)  plotted  versus  7*^  t a.(l  -+  B)/S. 


45 


144.2  kPa  (20.9  psi) 

Shot  l4  - 2 -Dimensional 
E(;uatlons  1 and  8 

r*  = t a^(l  + E)/3 


Figure  23  Average  front  surface  pressure  ratio  of  the  tvo-diaeneio&al 
shock  tube  target  for  an  incident  shock  overpreasure  of 
144,2  kBa  (20.9  P»l)  plotted  versus  •*  1 1 (l  ♦ B)/S. 


Figure  25  Avc-rage  surface  pressure  ratio  of  tbe  three-dttaenstoaal 
shock  tube  target  for  an  lucident  shock  overpressure  of 
73»7  kBa  (lO.T  pai)  plotted  versus  )***  « t a U + B)/^. 


1.1 


I 

'£ 


Figtm  26  Average  front  eurface  presaure  ratio  of  the  three-dijeienaioiial 
atoock  t«l»e  target  for  an  incident  eho&k  overpi*e88ure  of 
lU.l  {l6.l  pai)  plotted  veraua  >^  t a (1  + H)/8. 


Ilgur«  27  Average  front  aurface  presaure  ratio  of  the  three -ditaanaional 
ihock  tube  target  for  an  incident  ahoek  overpreaaure  of 
ikT.l  kfa  (21.3  psi)  plotted  versus  “ t a ( 1 B)/S. 
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One  requirement  for  use  with  the  overturning  code  is  that  the  predic- 
tion equation  used  should  provide  reasonable  predictions  over  the  range 
of  incident  shock  overpressures  from  13  to  345  kPa  (2  to  50  psi) . 

Figures  28  and  29  show  predictions  made  using  Equation  1 with  functions 
as  defined  by  Equation  4 and  Equation  8 for  a target  with  R = 0.5, 

S = 2 metres,  = 101.325  kPa  (14.969  psi)  and  = 15“C.  It  was 

assumed  that  no  decay  occurred  behind  the  shock  front . For  both  over- 
pressure extremes  the  fit  versus  r*  (Equation  8)  is  higher  initially  and 
lower  later  than  the  fit  versus  x (Equation  4).  No  data  were  available 
to  provide  a basis  for  selecting  one  fit  in  preference  to  another.  The 
fit  provided  by  Equations  1 and  4 was  selected  as  the  preferred  fit 
because  the  pressures  predicted  are  lower  at  early  times  than  the  fit 
provided  by  Equations  1 and  8.  Additional  data  are  required  to  determine 
which  of  the  two  fits  is  best,  and  whether  fitting  against  x rather  than 
agaiiist  will  provide  the  simplest  representation  over  the  widest  range 
of  incident  shock  overpressures. 


III.  REAR  SURFACE  LOADING 


In  Reference  2,  Taylor  presents  average  loading  data  for  rear  sur- 
faces for  the  configurations  shown  in  Figures  2C  and  2D.  Figures  30  and 
31  show  these  data.  The  average  pressure  ratio  for  the  rear  surface  is 
shown  plotted  versus  x-,  where  x„  = t a /S.  Note  that  this  time  normal - 

l'  p 0 

ization  does  not  include  the  parameter  R.  The  two-dimensional  configura- 
tion (H  s 0)  is  configuration  D,  and  the  throe-dimensional  configuration 
(R  ® 1)  is  configuration  C.  The  curves  were  obtained  by  averaging  the 
output  of  the  throe  gages  shown  in  Figure  I according  to  the  relation 


P 

average 


(2A  + C ♦ D)/4 


Although  this  relation  was  shown  to  provide  a good  average  for  the 
three-dimensional  configuration  at  .34.47  kPa  (5  psi)  in  Reference  .5, 
the  placement  of  gages  and  this  averaging  procedure  seems  less  suitable 
for  the  two-dimensional  configuration.  However,  there  wore  no  hydrocodo 
data  at  h.and  to  use  to  evaluate  this  averaging  procedure  for  the  two- 
dimensional  configuration. 

The  curves  from  Reference  2 shown  in  Figures  30  and  31  do  not  indi- 
cate the  initial  time  delay  which  corresponds  to  the  time  between  shock 
arrival  at  the  edge  of  the  rear  face  when  loading  of  the  rear  surface 
begins,  and  wave  arrival  at  the  gage  nearest  an  edge.  This  time  differ- 
ence was  not  measured  in  the  program  of  Reference  2.  However,  the  ini- 
tial wave  crossing  the  rear  surface  must  be  a weak  shock  wave  as  shown 
by  the  low  pressure  at  the  initial  rise  of  gage  records.  The  velocity 
of  this  initial  wave  must  be  very  near  the  sound  velocity  in  undisturbed 
air,  and  hence  this  sound  x^olocity  was  used  to  compute  the  initial  delay 
for  the  average  loading  records.  The  time  delays  expressed  in  terms  of 
Yg  were  0.488  for  the  two-ditsensionai  configuration  and  0.166  for  the 

three-dimensional  configuration. 

4$ 
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Equations  1 am  -t 
Eqtiations  1 and  8 

R » 0.5,  S = 2 metres 

= 101.325  kPa  (14.696  psi) 

=.  288.15  °K  (15  °C) 

P„  ==  13.78  kPa  (2.0  psi) 
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Figure  28  Predictions  of  froajt  eurfstce  pressure  ratios  using 
Equations  1 and  4 and  £<;Luationd  1 and  8 for  an 
incident  aUocH  overpressure  of  13«78  kBa  (2.0  pal) 


Figures  32  and  33  show  the  curves  from  Figu.es  30  and  3i  with  the  time 
delay  for  wave  arrival  included  so  that  fg  = 0 corresponds  to  the  initia- 
tion of  loading  on  the  rear  surface.  Values  were  estimated  for  the  aver- 
age maximum  pressure  ratio  achieved  on  the  rear  surface  at  the  conclusion 
of  the  initial  rise  of  the  record.  These  values  are  listed  as  P[y^^/Pg 

in  Table  1,  where  is  the  estimated  maximum  average  pressure  on  the 


RM 

rear  surface  and  is  the  incident  peak  overpressure, 
shown  plotted  versus  z » P /P  in  Figure  34. 

S A 


This  ratio  is 


A.  Fitting  an  Equation  for  Average  Rear  Surface  Pressure 


The  curves  in  Figures  32  and  33  for  experimental  data  were  normalized 
so  that  the  ordinate  value  of  1.0  corresponds  to  the  value  of 

was  found  that  the  following  relation  provided  what  seems  to  be  a good 
average  curve  for  both  the  two  and  three-dimensional  configurations: 


/-Ml  - (9) 

Figure  35  shows  the  norma’izod  curves  superimposed  without  regard  for 
individvuil  shot  identification  for  the  two  and  three-dimensional  configur- 
ations. The  dashed  curve  show*:  in  each  plot  is  the  plot  of  Equation  9. 


Figure  .36  shows  the  normalized  curves  for  the  two-dimensional  config- 
uration for  Shot  S (32.82  kPa  or  4.76  psi)  and  Shot  6 (139.4  kPa  or  20.22 
psi)  and  for  Equation  9.  l-iguro  3?  shows  a similar  plot  for  Shot  2 
(32.%  kPa  or  4.78  psi)  and  Shot  10  (148.S  kPa  or  21.54  psi)  fot  the 
three-dimensional  configuration,  lltere  Is  an  indication  of  a variation 
in  shape  with  z,  the  incident  shock  pressure  ratio.  In  Figure  37  the 
curve  for  the  lower  pressure  seems  to  rise  faster  than  the  curve  for  the 
higher  pressure.  The  curves  shown  in  Figure  36  aLso  suggest  a variation 
with  z.  in  the  initial  part  of  the  record  the  variation  is  the  opposite 
of  that  shown  in  Flgut'c  57,  while  for  tg  > 1.75  the  curve  for  the  higher 

£ '^oes  seem  to  be  lower  than  that  for  the  lower  z. 

Another  stiuly^  ^hows  a variation  wit.i  incident  overpressure  level  for 
a configuration  corresponding  to  R « I or  the  three-dimensional  configura- 
tion, Considering  this  fact,  and  that  the  more  usual  structure  of  inter- 
est more  nearly  corresponds  to  the  three-diaiensioiial  rather  then  the  two- 
dimensional  configuration,  an  equation  was  fitted  considering  the  2 
dependence  as  shown  by  the  three-dimensional  configuration.  Because  of 
the  uncertainty  in  the  accuracy  of  the  average  pressure  for  the  two-dimen- 
sional configuration  and  the  mixed  variation  with  z shewn,  no  attempt  was 
made  to  fit  a variation  with  z for  the  tifo-dimensionai  configuration* 
Instead,  the  same  variation  as  found  for  the  three-dimensional  configura- 
tion was  asauoed  to  apply  for  both  coses,  and  for  all  values  of  R betweim 
0 and  1. 
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vavQ  arrival  included.  Ploto  of  a fitted  equation  and  a 
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Figure  33  Average  preaeure  ratios  on  the  rear  face  of  a tltt*ee- 

dliaenslonal  ahoclt  tubs  target  plotted-  tneraus  “ t a^/S, 

vltb  delay  for  wave  arrlral  Included.  Plots  of  a fitted 
equation  and  a staudard  prediotlpn  technique  are  shown. 
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Figure  Average  maxlmUKi  pressure  ratio  on  rear  surface  of  shock 
tube  target  after  initial  pressure  rise  veroua  z 
Flots  of  a fitted  equation  are  shown.  ^ 
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Figure  35  Pr&sswe  on  the  rear  ourface  of  shock  tube  target 

ticumlised  to  average  maxiaum  pressure  after  Initial 
pressure  rise  versus  ^ “ t a /S,  Flota  of  a fitted 
equation  arc  showt^. 


I Figui’e  36  Comparison  of  pressure  ratios  on  rear  surface  of  a two- 

l dimensional  shock  tube  target  for  maximum  and  minimum 

I incident  shock  overpress\ires  used  in  an  experimental 

study*^.  The  plot  of  a fitted  equation  is  shown. 
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Figure  37  Coaj>arii?on  of  pressure  ratios  on  rear  surface  of  a three- 
dimensional  shock  tube  target  for  maximum  and  minimum 
ineidfiut  shock  overpreasitres  used  In  epn  ekperimentaX 
etudy*.  The  plot  of  a fitted  equation  is  shown. 


The  equation  fitted  was 


_L_  = 1 . * O.ldtg^)  (10) 

’’rm 

Figures  38  and  39  show  the  normalized  records  with  the  fit  from 
Equation  10  plotted  as  the  dashed  curve. 

To  develop  a complete  prediction  equation  for  the  rear  surface  pres- 
sure, it  was  necessary  to  find  equations  for  the  maximum  rear  surface 
pressure  used  as  normalizing  factors.  These  data  are  shown  in  Figure  34. 
A difficulty  added  in  the  fitting  process  was  the  requirement  for  finding 
relations  which  would  behave  in  a reasonable  fashion  when  used  outside 
the  range  of  fitting,  say  over  the  range  of  z from  0.1  to  3.4.  The  rela- 
tions found  were  as  follows: 

Two-dimensional,  R =>  o 


-0.73Z 


0.5  ♦ 1 .40 


-0.83Z 


Three-dimensional,  R » 1 


ill) 


''rH  ^ .0.32,“-5  ‘ 

s 


-0.83t 


(12) 


The  equations  where  combined  to  predict  values  for  the  variation  of 
R between  0 and  1 as  follows; 


’^RM  ^ -(0.73  - 0.4111)2^*®  * 


-0.852 


(13) 


The  curves  shown  in  Figure  54  are  plots  of  this  equation  with  R «•  0. 

The  final  prediction  equation  for  average  pressure  on  the  rear  surface 
is  the  coiabination  of  Equations  10  and  15.  and  is  as  follows: 


Jt 

P. 


73-  0.41R)2 


0.S+  l.4e 


•o»a5s 


1 -e 


-0.940 


*80  * 


(M) 


Figure  3S  Average  nomdlssed  treasure  ratio  on  tAie  rear  surface 
of  a three^>dijaen3loi^  shock  tube  target  ccsn^ed  vl^ 
a fitted  e<iviatica. 
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The  dashed  curves  in  Figures  32  and  33  were  calculated  using  this 
equation,  and  the  quality  of  the  fit  can  be  judged  from  these  plots.  The 
fit  is  best  for  the  three-dimensional  configuration,  where  R = 1. 

B.  Comparison  with  Other  Shock  Tube  Data 

In  Reference  6 rear  surface  average  pressures  are  shown  for  a config- 
uration for  R = 1.  Figure  40  shows  these  data  compared  with  the  predic- 
tions from  Equation  14.  The  multiple  curvatures  shown  by  the  experimental 
data  are  not  modeled  by  the  equation,  but  the  equation  does  seem  to  pro- 
vide a curve  which  fairs  reasonably  through  the  data. 

C.  Comparison  with  a Three-Dimensional  Hydrocode  Calculation 

Figure  41  presents  a comparison  of  the  predicted  curve  calculated 
using  Equation  14  with  results  from  a hydrocode  calculation  for  the 
average  pressure  on  the  rear  surface  of  an  S-280  Electrical  Equipment 
Shelter^.  The  calculation  was  made  for  an  incident  shock  front  over- 
pressure of  34.475  kPa  (S  psi) . It  was  assumed  that  no  decay  occurred 
behind  the  shock  front.  The  calculation  is  reported  in  detail  in 
Reference  7.  The  agreement  shown  in  Figure  41  seems  satisfactory. 

D.  Comparison  with  Standard  Prediction  Techniques 

In  Roferonco  3 the  average  pressure  rise  on  the  rear  surface  is 
approximated  by  a linear  increase  from  the  time  of  arrival  of  tho  shock 
front  at  tho  roar  surface  to  a maximum  value  given  by: 


^ i [1  ♦ (1  - (IS) 

s ^ 

The  time  for  tho  rise  is  specified  as  t • corresponds 

to  Tg  « 4. 

Predictions  using  this  technique  arc  shown  as  the  straight  iir.es  in 
Figures  32  and  S3,  with  the  X indicating  the  maximum  value,  live  large 
discrepancy  emphasized  by  Taylor^  between  the  straight-line  approximation 
and  the  shock  tube  data  is  evident. 

b-  Modification  to  Account  for  Decay  of  the  Incident  Blast  Wave 

When  the  overpressure  behind  the  incident  shock  front  decreases  with 
time  as  occurs  in  a blast  wave,  the  average  pressure  on  tho  rear  surface 
is  assumed  to  be  reduced  in  proportion  to  the  reduction  in  incident  pres- 
sure from  its  peak  value.  Ei|uation  14  remains  unchanged  except  that  P^(t) 

is  substituted  for  P,  the  peak  Incident  uverpressuro. 
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Figure  Uo  Coraparison  of  rear  surface  pressure  ratios  for  the  three- 
dimensional  shock  tulje  target  of  Reference  6 vlth  an 
equation  fitted  to  the  ehock  tube  data  of  Reference  2. 
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Modification  for  Code  Use 


The  suggested  form  of  the  average  front  surface  pressure  function 
for  incorporation  into  the  overturning  code  of  Reference  1 decays  to 
stagnation  overpressure  as  shown  in  Equation  4.  The  particular  over- 
turning code  terminates  the  calculation  of  diffraction  loading  when  the 
overturning  moment  due  to  drag  loading  becomes  equal  to  or  exceeds  that 
due  to  the  diffraction  moment.  Assuming  equal  diffraction  and  drag 
loading  areas  and  moment  arms,  the  moments  are  equal  at  late  times  when: 

stag  rear  d q 
where:  ^stag  " stagnation  overpressure 

^rear  * average  overpressure  on  rear  surface 

= drag  coefficient 
P^  = peak  dynamic  pressure. 

To  insure  that  the  diffraction  loading  calculation  terminates  in  a 
reasonable  time  with  an  appropriate  match  between  the  diffraction  and 
drag  loading,  the  asymptotic  value  of  P ^ was  adjusted  sc  that 


> - P = 0.97  C,  P . 

stag  rear  d q 


The  choice  of  the  factor  0.97  is  arbitrary,  but  its  use  is  expected 

to  produce  a transition  where  the  diffraction  loading  essentially  has 

been  completed.  Then  the  value  for  P is: 

^ roar 


P = P , 
roar  stag 


0.97  C.  P„. 

d q 


(lb) 


For  a Hast  wave  where  the  pt'ossuro  decays  with  time,  and 

P^{t)  are  substituted  for  a«d  respectively.  Then  the  predic- 

tion relation  appropriate  for  incorporation  into  an  overturning  code  ;is: 


.-0.94e‘°*^^*Tft(i  ♦ O.lbt  ^ 


i 


(i7) 


IV.  Dia'USSION 

fiqujjtions  for  predicting  average  pressure  on  the  front  and  rear 
surfaces  of  a rectangular  paraUclepiiiod  have  been  developed  by  fairing 
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curves  through  a limited  set  of  data  from  Reference  2.  The  equations 
fitted  for  the  three-dimensional  configuration  provide  reasonable  agree- 
ment with  shock  tube  data  from  References  2 and  6 in  the  overpressure 
range  from  22  to  138  kPa  (3  to  20  psi)  and  with  a hydrocode  calculation^ 
at  34.475  kPa  (5  psi).  The  validity  of  the  fits  for  the  two-dimensional 
configuration  is  subject  to  question  because  the  placement  and  weighting 
of  the  gage  records  in  Reference  2 does  not  seen;  as  suitable  for  the 
two-dimensional  configuration  as  for  the  three-dimensional  configuration. 
Thus  the  equations  fitted  for  the  two-dime.isional  configuration  should 
be  regarded  as  an  interim  fit  which  may  be  improved  if  additional  shock 
tube  data  and  hydrocode  calculations  become  available  and  are  utilized. 
Because  the  two-dimensional  results  are  incorporated  in  the  final  pre- 
diction equations  where  R varies  between  0 and  1,  the  confidence  in  the 
prediction  should  decrease  somewhat  as  R approaches  zero. 

;.t  was  found  that  the  available  shock  tube  data  for  average  front  sur- 
face overpressure  could  be  fitted  satisfactorily  against  either  the  dimen- 
sionless pai'ameter  t = t a^Cl  + R)/S  or  x*  * t (1  R)/S.  Addi.ional 

data  for  a wider  range  of  shock  overpressures  are  required  to  determine 
which  of  the  two  parameters  can  be  used  to  obtain  the  most  useful  fit. 

The  prediction  equations  provide  smooth  extrapolations  to  the  extreme 
values  of  overpressure  of  interest  (13  to  345  kPa,  or  2 to  50  psi),  but 
the  validity  of  the  extrapolation  at  these  extremes  is  not  established. 
Additional  data,  particularly  at  overpressures  above  20  psi,  are  required 
for  values  of  R between  0 and  1 to  provide  a basis  for  properly  accounting 
for  the  variation  of  predicted  average  pressure  with  R Tlie  modifications 
to  account  for  decay  of  pressure  behind  the  incident  shock  fixnt  should 
bo  evaluated  if  appropriate  data  become  available. 
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LIST  OF  SYMBOLS 


Symbol 


A 


a 

a 

0 


B 


= function  of  time,  incident  overpressure,  ambient  overpres- 
sure, and  R used  in  equations  fitted  to  average  pressure 
profiles. 


= parameter  used  in  a fitted  equation. 


= ambient  sound  velocity  = 331.32 


v; 


Ta 


273.16 


m/s. 


= velocity  of  the  rarefaction  wave  behind  the  reflected  shock 


wave  = a 


(Sz  7)(2z  7; 

7(7  + 6zj 


= function  of  time,  incident  overpressure,  ambient  overpres- 
sure, and  R used  in  equations  fitted  to  average  pressure 
profiles. 


= drag  coefficient  for  a blast  target. 


CONF  = symbol  meaning  configuration 

P = average  overpressure  on  front  surface  or  roar  surface. 

P(t)  = average  overpressure  versus  time  for  a front  or  rear  surface, 

^average  ’ pressure  on  front  surface. 


= ambient  atmospheric  pressure. 


overpressure  approached  as  a limit  on  the  front  surface. 


P B peak  dynamic  pressure. 


Pq(t) 


o dyiiaroic  pressure  versus  time  for  a blast  wove. 


Pn  o Peak  reflected  overprossuia  c 2 P,  (4t  ♦ 7)/(2  ♦?). 

K S 


P « average  pressure  on  roar  surface  approached  as  a limit  as 

time  increases . 

Pfy.  o estimated  maximum  average  overpressure  on  rear  surface  after 
initial  pressure  rise. 
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LIST  OF  SYMBOLS  (Continued) 


Symbol 

P 

s 

P,.(t) 

p 

stag 


= peak  incident  shock  overpressure. 

= overpressure  versus  time  for  a blast  wave. 

= stagnation  overpressure  on  the  front  surface;  which  for  a 
non-decaying  shock  wave  and  a Mach  number  less  than  one  is: 

2 \7/2 


stag 


(P  ♦ 
s 


* 7(2  * 1)(Z  W)  ) 


- P. 


stag 


(t)  >=  stagnation  overpressure  at  time  t on  front  surface,  which 


t* 


<> 

B 

y 

T 


for  Mach  number  less  than  one  is: 


P * (t)  = (P  (tl  + PJ 

stag  ' s A 


1 * (2/7) 


1 


7/2 


- P 


A 


= ratio  of  the  edges  of  the  minimiua  I'octangular  area  which 
roust  be  examined  to  define  loading  over  the  entire  surface 
area.  The  ratio  is  taken  so  that  R ranges  from  0 for  tlie 
two-dimensional  configuration  to  I for  the  throe*dimonsional 
configuration. 

“ roinimum  dlstaj\ce  of  travel  from  an  edge  of  a front  or  rear 
surface  for  a wuvo  to  traverse  tho  surface. 

« time, 

e time  for  decay  of  front  surface  average  overpressure  in  a 
standard  prediction  technique^,  t^  “ 3 ^^^5* 

« ambient  atmo.sphoric  toraperaturo,  degrees  Kelvin. 

“ ratio  of  peak  incident  shock  overpressure  to 

arobient  atjaospheric  pressu2’e. 

» parameter  used  in  e fitted  equation. 

» parameter  used  in  a fitted  equation. 

« ratio  of  specific  heats  y ■ 1.4  for  an  ideal  gas. 

« t a^d  ♦ R)/S. 
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LIST  OF  SYMBOLS  (Continued) 


I 

y 


Symbol 

T*  = t a^(l  + R)/s. 

f Tb  = t a^/S. 

2DF  = s>'mbol  meaning  two-dimensional  configuration,  front  surface 
record . 

2DR  = symbol  meaning  two-dimensional  configuration,  rear  surface 
record . 

SDF  = symbol  meaning  three-dimensional  configuration,  front  sur- 
face record. 

3DR  o symbol  meaning  three-dimensional  configuration,  roar  surface 
record . 
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ATTN:  Doc  Control  for  Tech  Lib 
P.  0.  Box  969 
Livermore,  CA  94550 


2 TRW  Systems  Group 

ATTN:  Benjamin  Sussholtz 

Tech  Info  Ctr/S-I930 
One  Space  Park 
Redondo  Beach,  CA  90278 


1 Science  Applications,  Inc. 
ATTN:  Teclinical  Library 
P,  0.  Box  3507 
Albuquerque,  NM  87110 

2 Science  Applications,  Inc 
ATTN:  R.  Seebaugh 

John  Mansfield 
1651  Old  Meadow  Road 
McLean,  VA  22101 

1 Science  Applications,  Inc. 
8201  CapW'jll  Drive 
Oakland,  CA  94621 

2 Science  Applications,  Inc. 
ATTN : Technical  Library 

Michael  McKay 
P.  0.  Box  2351 
La  Jolla,  CA  92038 

2 Systems,  Science  ^ Software 
ATTN;  Donald  R.  Grine 
Technical  Library 
P.  0.  Box  1620 
La  Jolla,  CA  92037 

1 Terra  Tek,  Inc. 

ATTN:  Technical  Library 

420  Wakara  Way 

Salt  Lake  City,  UT  84108 

1 Totra  Tech,  Inc. 

ATTN:  Technical  Library 
630  North  Rosemead  Blvd. 
Pasadena,  CA  91107 


1 TRW  Systems  Group 
ATTN:  Greg  Hulcher 

San  Bernardino  Operations 

P.  0.  Box  1310 

San  Bernardino,  CA  92402 

2 Union  Carbide  Corporation 
Holifield  NationaT  Laboratory 
ATTN:  Doc  Control  for 

Tech  Lib 

Civil  Defense  Research 
Proj 

P.  0.  Box  X 

Oak  Ridge,  TN  37830 

1 Universal  Analytics,  Inc. 
ATTN:  E.  I.  Field 
7740  W.  Manchester  Blvd. 

Playa  del  Rey,  Ca  90291 

1 Woidlingor  Assoc.  Consulting 
Engineers 

ATTN:  M.  L.  Baron 
110  East  59th  Street 
Now  York.  NY  10022 

1 Wostinghouse  Electric  Company 
Marine  Division 

ATTN;  W.  A.  Votz 
liendy  Avenue 
Sunnyvale,  CA  94008 

2 Battolle  Memorial  Institute 
ATTN:  Technical  Library 

R.  W.  Klingosmith 
505  King  Avenue 
Columbus,  OH  43201 
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1  California  Institute  of 
Technology 
ATTN:  T.  J.  Ahrens 
1201  E.  California  Blvd. 
Pasadena,  CA  91109 

1 COSMIC 

ATTN;  L.  C.  Gadol 
112  Barrow  Hall 
University  of  Georgia 
Athena,  GA  30602 


No.  of 

Copies  Organization 

2  Southwest  Research  Institute 
ATTN:  Dr.  W.  E.  Baker 
A.  B.  Wenzel 
8500  Culebra  Road 
San  Antonio,  TX  78206 

2 Stanford  Research  institute 
ATTN:  Dr.  G.  R.  Abrahamson 
Carl  Peterson 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 


2 Denver  Research  Institute 
University  of  Denver 
ATTN:  Mr.  J.  Wisotski 

Technical  Library 
P.  0.  Box  10127 
Denver,  CO  80210 

2 IIT  Research  Institute 
ATTN:  Milton  R.  Johnson 

Technical  Library 
10  West  3Sth  Street 
Chicago,  IL  o0616 

3 Lovelace  Foundation  for 

Medical  Education 
ATTN:  Asst.  Dir.  of  Reseax’ch/ 
Robert  K.  Jones 
Technical  Library 
Dr.  D.  Richmond 
5200  Gibson  Blvd,  SE 
Albuquerque,  NM  87108 

1 Massachusetts  Institute  of 
Technology 

Aeroelastic  and  Structures 
Research  Laboratory 
AITN:  Dr.  E.  A.  Witmor 
Cambridge,  MA  02139 


1 University  of  Illinois 
Consulting  Engineering  Services 
ATTN:  Nathan  M.  Newmark 

1211  Civil  Engineering  Building 
Urbana,  IL  61801 

2 The  University  of  New  Mexico 
The  Eric  H.  Wang  Civil 

Engineering  Research  Facility 
ATTN:  Larry  bickle 
Neal  Baum 
University  Station 
Box  188 

Albuquerque,  NM  87131 

2 Washington  State  University 
Administration  Office 
ATTN:  Arthur  Miles  Hohorf 
George  Duval 
Pullnan,  WA  99163 

Aberdeen  Proving  Ground 

Marine  Corps  Ln  Ofc 
Dir,  USAMSAA 
ATTN:  Dr.  J.  Sperrazza 
Mr.  R.  Norman.  GWD 
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